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Abstract 

Acid rain (AR) is a serious environmental issue inducing harmful impacts on plant growth and development. It has been 
reported that Liquidambar formosana, considered as an AR-sensitive tree species, was largely injured by AR, compared with 
Schima superba, an AR-tolerant tree species. To clarify the different responses of these two species to AR, a comparative 
proteomic analysis was conducted in this study. More than 1000 protein spots were reproducibly detected on two- 
dimensional electrophoresis gels. Among them, 74 protein spots from L. formosana gels and 34 protein spots from S. 
superba gels showed significant changes in their abundances under AR stress. In both L formosana and S. superba, the 
majority proteins with more than 2 fold changes were involved in photosynthesis and energy production, followed by 
material metabolism, stress and defense, transcription, post-translational and modification, and signal transduction. In 
contrast with L. formosana, no hormone response-related protein was found in S. superba. Moreover, the changes of 
proteins involved in photosynthesis, starch synthesis, and translation were distinctly different between L. formosana and S. 
superba. Protein expression analysis of three proteins (ribulose-1,5-bisphosphate carboxylase/oxygenase large subunit, 
ascorbate peroxidase and glutathione-S-transferase) by Western blot was well correlated with the results of proteomics. In 
conclusion, our study provides new insights into AR stress responses in woody plants and clarifies the differences in 
strategies to cope with AR between L. formosana and S. superba. 
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Introduction 

Acid rain (AR) emerged as a serious environmental issue as a 
consequence of the increasing industrial activities throughout the 
world [1]. Forty percent of the territory in China is seriously 
affected by AR since the late 1970s, especially in southern China 
[2] . The harmful impacts of AR on plants are observed in a wide 
array of biological processes. AR decreases seed germination [3], 
strips the protective wax from leaves [4], induces visible injury 
symptoms [1], disturbs plant nitrogen metabolism [5]. AR also 
decreases chlorophyll content and photosynthetic efficiency [3,6], 
increases reactive oxygen species (ROS) production [7], acceler- 
ates the leaching of nutrients from plant foliage [8] , which further 
inhibits tree radial growth, vertical growth and total tree biomass 
[4,9]- ^ 

Liquidambar formosana and Schima superba are both dominant 
broad-leaf tree species and are distributed over large surface areas 



in the forest of southern China [10]. Some field observations and 
laboratory experiments reported that, when compared with S. 
superba, L. formosana was largely injured by AR during the past 
decade, which had negative impacts on forest ecosystem [9]. Our 
previous study also found that AR more easily affected some 
physiological parameters in L. formosana seedlings than S. superba's 
seedlings, e.g., seed germination, seedling growth, photosynthesis, 
antioxidant system, etc. [3,1 1]. Thus, L. formosana is considered as 
an AR-sensitive species, while S. superba is an AR-tolerant species. 
Although the differential responses of L. formosana and S. superba to 
AR have been analyzed at the morphological and physiological 
level, a comprehensive elucidation of the molecular mechanisms 
underlying the different strategies to cope with AR between two 
tree species is still needed. 

Proteomics is a powerful tool for providing new insights into 
complete proteomes at the organ, tissue and cell levels [12]. A 
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number of proteomic analyses help us understand the molecular 
mechanisms of plants in responses to various environmental 
stresses including salinity [13], cold [14], heavy metal [15], etc. In 
our previous work, a wide array of proteins related to AR- 
resistance has been identified by comparative proteomic analysis 
in a model plant, Arabidopsis thaliana [5,16]. However, little 
information is available in proteome analysis for tree species 
subjected to AR stress, and a comparison between AR-sensitive 
and AR-resistant broad-leaf tree species has not been fully 
conducted at the proteome level. 

In the present study, we initiated a comparative proteomic study 
to systematically investigate the changes in protein profile in two 
broad-leaf tree species, L.formosana and S. superba, that are different 
sensitive to AR tolerance, when submitted to simulated AR 
(pH 3.0) treatment for one month. Based on two-dimensional 
electrophoresis (2-DE) and mass spectrometry (MS) analysis, a 
comprehensive inventory of proteins regulated by AR was 
established in the two tree species. The overall objectives of this 
study are (1) to provide valuable insights into AR stress responses 
in woody plants; (2) to clarify the differences in strategies to cope 
with AR between L. formosana and S. superba. 

Materials and Methods 

Plant materials and treatments 

Seeds of L. formosana and S. superba were purchased from Tree 
Seed Centre of Shuyang County in Jiangsu Province, China. The 
seeds have been mixed together when they were collected from 
independent individuals and families. Seeds were surface-sterilized 
with 0.5% hypochlorite for 30 min, then washed thoroughly 
with distilled water. Then the seeds were germinated in a soil/ 
vermiculite (1:1) mixture in an environmentally controlled growth 
chamber. For each species, three weeks old healthy seedlings with 
similar size were randomly transplanted into individual pots, each 
with a dimension of 24 cm (open top) x 1 3 cm (height) x 1 5 cm 
(flat bottom), and filled with soil/vermiculite (1:1) mixture. Fifteen 
seedlings were planted in one pot. All seedlings were cultivated in 
the same controlled growth chamber with a daily temperature 
regime of 28/25°C (day/night), relative humidity of 60-70% and 
a 12-h photoperiod at 210 |imol m -2 s _1 photosynthetically active 
radiation (PAR). Three months later, the seedlings were divided 
into control group (GK) and simulated AR treatment group (AR). 
Each group had at least three replicates. The control group and 
simulated AR treatment group were sprayed once per day with the 
control (pH 5.6) solution and AR (pH 3.0) solution, respectively. 
The ion compositions of the control solution was adapted from 
Fan and Wang [4], AR solution was made from control solution 
and the pH was adjusted by adding a mixture of H 2 S0 4 and 
HN0 3 in the ratio of 5:1. The final concentration of H2SO4 and 
HN0 3 were 0.45 and 0.09 mM, respectively, which represents 
the average ion composition of rainfall in southern China [4]. 
After one month of treatment, a portion of fresh leaves was used 
for measuring some physiological parameters such as necrosis 
percentage, chlorophyll content, net photosynthetic rate, H 2 0 2 
content and so on. The remaining leaves were stored at — 80°C for 
proteomic and Western blot analysis. 

Necrosis percentage and chlorophyll content 
measurements 

At least 30 fully expanded leaves of each species were randomly 
selected from control and AR treatment groups and photographed 
using a digital camera. The necrosis percentage was calculated as 
described previously [16]. Chlorophyll was extracted from 
approximately 0.1 g fresh leaf slices directly into 10 ml ice-cold 



acetone (80%, v/v). The chlorophyll content (mg g 1 fresh weight 
(FW)) was determined as described previously [1 7] . 

Net photosynthetic rate and chlorophyll fluorescence 
measurements 

Three seedlings per species were randomly chosen from 
different pots and at least two fully emerged leaves per plant 
were selected for net photosynthetic rate (Pn) and chlorophyll 
fluorescence measurements. Pn was performed with a portable 
photosynthesis system (LI-6400, Li-Cor Inc., Lincoln, Nebraska, 
USA), as described previously [18]. According to the method of 
Liu et al [16], leaf chlorophyll fluorescence was measured using 
a pluse-amplitude-modulation fluorometer (PAM-2100, Heinz 
Walz, Effeltrich, Germany). 

Measurement of proline, malondialdehyde (MDA) and 
ROS production 

Proline content was measured according to the method of 
Jiang et al [19]. The level of lipid peroxidation was measured by 
estimating MDA content using thiobarbituric acid (TBA) reaction 
[20]. Superoxide radical (0 2 '~) and H 2 0 2 content was measured 
following the method of Chen et al [1 1]. 

Protein extraction 

Protein extraction was performed using phenol-based protocol 
described by Liu et al [5], with slight modifications. Briefly, frozen 
leaves (1.0 g) were ground with a mortar and pesde with liquid 
nitrogen, the ground powder was homogenized in pre-cooled 
extraction buffer (20 mM Tris-HCl, pH 7.5, 250 mM sucrose, 
10 mM ethylene diamine tetraacetic acid (EDTA), 1 mM 
phenylmethyl-sulfonyl fluoride, 1 mM 1 ,4-dithiothreitol (DTT) 
and 1% Triton X-100) on ice. Then an equal volume of ice-cold 
Tris-HCl (pH 7.5) saturated phenol was added and the mixture 
was centrifuged (15,000 4°C) for 15 min. The phenol phase 
was collected and proteins were precipitated with ammonium 
acetate in methanol for 10 h at — 20°C. After centrifugation, the 
supernatant was discarded and the pellet was washed for three 
times using cold acetone containing 10 mM DTT. The washed 
protein pellets were dissolved in a lysis buffer (8 M urea, 2 M 
thiourea, 4% (w/v) 3-[(3-cholamidopropyl)dimethylammonio]- 
1-propane sulfonate (CHAPS), 1% (w/v) DTT and 0.5% (w/v) 
IPG buffer pH 4-7) at room temperature. The total protein 
concentration of the lysates was determined using a Bio-Rad 
protein assay reagent (Bio-Rad, Hercules, CA, USA). 

Two-dimensional electrophoresis, image and data 
analysis 

Two-dimensional electrophoresis was conducted according to 
the methods of Bai et al [13] and Hu et al [21]. Isoelectric 
focusing (IEF) was done using an Ettan IPGphor system (GE 
Healthcare) PROTEAN electrophoresis system and immobilized 
IPG dry gel strips with a linear pH range (18 cm long, pH 4-7 
linear) (GE Healthcare Amersham Bioscience, Little Chalfont, 
UK). Protein samples (800 ixg) were loaded during the rehydration 
step at room temperature for 12 h. IEF was performed at 300 
voltage (V), 500 V and 1,000 V for 1 h, a gradient to 8,000 V over 
4 h, and kept at 8,000 V for a total of 80,000 volt-hours (Vh) at 
20°C. Subsequently, focused strips were equilibrated in equilibra- 
tion buffer as described by Yang et al [12]. For the second 
dimension, proteins were separated on 15% SDS polyacrylamide 
gels. Proteins spots were detected by staining the gels with 
Coomassie Brilliant Blue R-250. The 2-DE gels were scanned with 
a scanner (Uniscan M3600, China) and the gel images were 



PLOS ONE I www.plosone.org 



2 



July 2014 | Volume 9 | Issue 7 | e102532 



Proteomics Reveal AR Tolerance Metabolisms in Two Trees 



analyzed with PDQuest software (Version 8.01, Bio-Racl, Hercu- 
les, CA), on the basis of their relative volume. Only those protein 
spots with significant (more than 2-fold change) and reproducible 
changes in three replicates were selected for next MS analysis. 

In-gel digestion, protein identification and classification 
analysis 

The protein spots, which were differentially displayed under 
AR treatment, were excised from the preparative 2-D gels and 
digested by trypsin. After digestion, the peptide solution was 
collected and peptide mass fingerprint (PMF) was acquired using 
matrix- assisted laser desorption/ionization-time-of-flight mass 
spectrometry (MALDI-TOF MS) analysis (ReFlexTM III, Bruker, 
Bremen, Germany) as described previously [16]. The PMF spectra 
were used in online searches combined with the Mascot program 
search engine (http://www.matrixscience.com) and National 
Center for Biotechnology Information (NCBI) protein database 
(http://www.ncbi.nlm.nih.gov) (NCBInr, 17751536 entries, 
downloaded on April 17, 2012). PMF search parameters were 
set up as described previously [22]. Proteins with a MOWSE score 
>73 were considered as positive identifications. The identified 
proteins were searched with against the UniProt (http:/ /www. 
uniprot.org) and/or NCBI protein database (http:/ /www.ncbi. 
nlm.nih.gov) for updated annotation and homologous proteins 
identification. Afterwards, the successfully identified proteins 
were further classified using Functional Catalogue software 
(http:/ / mips.gsf.de/ projects/ funcat). 

Western blot analysis 

Western blot analysis was performed as described previously 
[11]. Total proteins (40 ug) extracted from L. formosana. and S. 
superba leaves were separated by 12% w/v standard sodium 
dodecyl sulfate polyacrylamide gel electrophoresis and then 
electrophoretically blotted to polyvinylidene difluoride membrane 
for 50 min. The membranes were blocked over-night with 
Western Blocking Buffer (TIANGEN, China). Protein blots 
were probed with primary antibodies of Rubisco large subunit 
(RuBISCO LSU) (AS03037-200, Agrisera, Sweden), ascorbate 
peroxidase (APX) (AS08368, Agrisera, Sweden) and glutathione- 
S-transferase (GST) (AS09479, Agrisera, Sweden), at dilution of 
1:5000, 1:2000 and 1:1000, respectively, for 4 h at room 
temperature with agitation. Next, the membranes were washed 
in phosphate buffered saline with Tween-20 solution (PBST) 
solution (50 mM Tris-HCl, pH 8.0, 150 mil NaCl, 0.05% Tween 
20, v/v) three times and incubated with anti-rabbit IgG 
horseradish peroxidase (HRP) conjugated to alkaline phosphatase 
(Abeam, U.K., 1:5000 dilution) for 1 h at room temperature to 
detect primary antibodies. /?-actin (1 :5000, Santa Cruz, California, 
USA) was used as an internal control. After several washes with 
PBST solution, membranes were incubated in an enhanced 
chemiluminescence (ECL) substrate detection solution (TIAN- 
GEN, China) according to the manufacturer's instructions. Images 
of protein blots were obtained using a CCD imager (FluorSMax 
Bio-Rad, Hercules, CA, USA). The optical density values of the 
protein signals were quantified using the Quantity One software 
(Bio-Rad, Hercules, CA, USA). 

Statistical analysis 

For physiological measurements, at lease four independent 
repetitions were used. Values in figures and tables were expressed 
as means ± se. The statistical significance of the data was analyzed 
using a univariate analysis of variance (One-way ANOVA, 
Duncan's multiple range test, /><0.05) with the SPSS 20.0 package 



(SPSS, Chicago, Illionis USA). For proteomic experiment, protein 
samples for 2-DE gel image analysis were extracted from three 
independent seedlings grown in three different pots in the same 
growth chamber. Thus, for each species, three independent 
biological replicates were performed in 2-DE gel image analysis. 
Statistic analysis for protein spot on 2-DE gels was performed 
using Student's t-test (/><0.05) provided by PDQuest software as 
mentioned earlier. 

Results 

Morphological and physiological responses of L 
formosana and S. superba to AR 

As shown in Fig. 1A and B, remarkable yellowing symptom and 
significant necrosis emerge in L. formosana leaves after one-month 
exposure to AR. There was an important decrease in total 
chlorophyll content in AR- treated L. formosana leaves, however, no 
statistically significant change was found in S. superba (Fig. 1C). Pn 
and Fv/Fm in AR-treated L. formosana seedlings were remarkably 
inhibited, whereas AR slighdy decreased Pn and Fv/Fm in S. 
superba (Fig. ID and E). After AR treatment, proline content in 
L. formosana and S. superba increased by 76.0% and 19.7%, 
respectively, compared with the control (Fig. 2A). MDA contents 
in AR-treated L. formosana and S. superba increased by 89.4% and 
44.8%, respectively (Fig. 2B). As shown in Fig. 2C and D, the 
levels of H 2 O z and O z "~ in both L. formosana and S. superba were 
significandy stimulated by AR. In particular, compared with the 
control, H 2 0 2 and 0 2 " content increased by 83.3% and 67.8% in 
L. formosana, and by 38.4% and 44.7% in S. superba, respectively 
(Fig. 2C and D). 

Identification of AR-responsive proteins 

To investigate the differentially expressed proteins in L. 
formosana and S. superba exposed to AR treatment, comparative 
proteomic analysis was performed on Coomassie-stained 2-DE 
maps shown in Fig. 3. Over 1000 protein spots reproducibly 
separated and matched between control and AR gels, 74 protein 
spots in L. formosana, and only 34 protein spots in S. superba had at 
least a 2-fold greater abundance in either AR or control (Fig. 3). 
The identified proteins in L. formosana and S. superba by MALDI- 
TOF MS analysis are presented in Tables 1 and 2, respectively. 
Because the complete annotated sequences of L. formosana and S. 
superba genomes are not yet available, all identified proteins were 
functionally classified by UniProt and NCBI databases according 
to their homology with other proteins. Functional annotations in 
databases existed for the majority of the protein spots, while 12 
proteins (spots L63-L74) in L. formosana and 6 proteins (spots S29- 
S34) in S. superba were annotated as predicted or unknown proteins 
(Tables 1 and 2). 

Of the 74 protein spots identified in L. formosana, the abundances 
of 53 proteins were increased and those of 21 proteins were 
decreased in response to AR (Table 1). In S. superba, 21 proteins 
were increased in their abundances and 1 1 proteins were 
decreased under AR stress (Table 2). Among these affected 
proteins, phosphoglycerate kinase (PGK, spot L32, L8, S6) and 
ATP synthase CF1 alpha subunit (spot L22, SI 1), were decreased 
in their abundances in both L. formosana and S. superba (Tables 1 
and 2). Remarkably, after AR treatment, abundance of maturase 
K was increased in L. formosana, but decreased in S. superba 
(Tables 1 and 2). Further analysis on the results revealed that some 
proteins were represented by more than one spot. These proteins 
included phosphoglycerate kinase (PGK, spot L32, L8), 
ATP synthase beta subunits (spot L22, L26) and maturase K 
(spot L57-59) in L. formosana and ribulose-l,5-bisphosphate 
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Figure 1. Effects of one-month AR on morphology and photosynthesis of L. formosana and S. superba. The pH of AR solution was 
adjusted to 3.0 by adding a mixture of H 2 S0 4 and HN0 3 in the ratio of 5:1. The final concentration of H 2 S0 4 and HN0 3 were 0.45 and 0.09 mM, 
respectively. (A) Leaf injury phenotype. (B) Leaf necrosis percentage. (C) Total chlorophyll content. (D) Net photosynthetic rate (Pn). (E) Quantum 
efficiency of open PSII centers in a dark-adapted state (Fv/Fm). Columns labeled with different letters indicate significant differences at p<0.05. 
doi:10.1371/journal.pone.0102532.g001 



carboxylase/oxygenase (Rubisco) large subunit (spot S9, SI 3, SI 4) 
in S. superba. The multiple spots might represent isoforms or 
different post- translation modification of individual proteins [13]. 

Functional classification of AR-responsive proteins 

In order to obtain annotation of AR-responsive protein, 
all identified proteins were further classified according to their 
biological function and cellular component categories in UniProt 



(http://www.uniprot.org) and/or NCBI protein database (http:// 
www.ncbi.nlm.nih.gov). AR-responsive proteins were found to be 
involved in a wide range of biological processes. After AR 
treatment, with the exception of photosynthesis and energy 
production related proteins, the abundance of most proteins were 
decreased in both L. formosana and S. superba (Fig. 4). As shown 
in Fig. 5, a higher percentage of proteins were involved in 
photosynthesis and energy production, which accounted for 
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Figure 2. Changes in proline (A), MDA (B), H 2 0 2 (C) and 0 2 ' (D) content after AR treatment. Columns labeled with different letters 

indicate significant differences at p<0.05. 

doi:10.1371/journal.pone.0102532.g002 



24.3% and 29.4% of the total proteins in L. formosana and S. 
superba, respectively. The following groups were proteins involved 
in material metabolism, stress and defense, transcription, post- 
translational and modification, and signal transduction. In 
opposition to L. formosana, no protein related to hormone response 
was found in S. superba. 

Protein expression analysis by Western blot 

Our above proteomic results revealed that the abundance of 
Rubisco was decreased (spot L31, L34), while APX (spot L40) and 
GST (spot L41) were increased in L. formosana under AR treatment 
(Table 1). As shown in Fig. 6A and B, compared with the control, 
the protein expression of Rubisco large subunit analyzed 
by Western blot was significantly decreased in AR-treated L. 
formosana. In contrast, the protein expression level of APX and 
GST increased 1.3-fold and 1.6-fold, respectively, compared to 
control, in L. formosana (Fig. 6C and D). No significant change in 
the expression of Rubisco large subunit, APX and GST was 
observed in S. superba (Fig. 6). 

Discussion 

AR has negative effects on plant growth and development [1,6]. 
Neves et al [23] found that simulated AR (pH 3.1) caused 
chlorosis and necrosis in leaves and led to chlorophyll loss and 
photosynthetic depression in Eugenia uniflora. Moreover, chloro- 
phyll content and photosynthesis in L. formosana were also 
remarkably suppressed by AR treatment in this study (Fig. 1). 
Compared with S. superba, reductions on chlorophyll content and 



photosynthetic ability by AR were more obvious in L. formosana 
(Fig. 1), which is consistent with the results of previous studies 
[9,11]. Besides chlorophyll content and photosynthesis, proline 
content, MDA content and ROS (e.g. H 2 0 2 and 0 2 " ) production 
are commonly used as biochemical markers to monitor the 
damage level in plants under environmental stress [24]. In this 
study, AR increased proline content, MDA content and ROS 
production in both L. formosana and S. superba (Fig. 2), which were 
consistent with the results obtained by Chen et al [11]. However, 
the increase in these physiological changes caused by AR was less 
pronounced in S. superba than those in L. formosana (Fig. 1 and 2), 
suggesting that S. superba, a tolerant species, had less cell damage 
than L. formosana, a sensitive species. 

To further reveal the different strategies to cope with AR 
between the two species, 74 protein spots in L. formosana and 34 
protein spots in S. superba caused by AR were identified by 
proteomic analysis in this study. Interestingly, similar results were 
also reported in previous studies, which reported more changes in 
protein abundance in sensitive species, Arabidopsis thaliana, than in 
tolerant species, Thellungiella halophila, under salt stress [25]. Since 
L. formosana had much higher changes in its protein profile, our 
results proved that this species is more sensitive to AR than S. 
superba. 

Photosynthesis and energy production-related proteins 

Photosynthesis is an essential metabolic process of plants and is 
vulnerable to environmental stress. It is well known that AR can 
remarkably reduce photosynthesis [3,23]. In this study, two light 
reaction-related proteins, including chlorophyllide a oxygenase 
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Figure 3. 2D gel analysis of proteins extracted from L. formosana and S. superba leaves. The numbers assigned to the proteins spots 
correspond to those listed in Tables 1 and 2. (A) Representative 2-DE gels of L formosana in which 74 spots showing at least 2-fold changes (p<0.05) 
under AR were identified by MALDI-TOF MS. (B) Close-up views of differentially expressed protein spots in L formosana (highlighted by arrows). (C) 
Representative 2-DE gels of S. superba in which 34 spots showing at least 2-fold changes (p<0.05) under AR were identified by MALDI-TOF MS. (D) 
Close-up views of differentially expressed protein spots in S. superba (highlighted by arrows). 
doi:10.1371/journal.pone.0102532.g003 



(CAO, spot L21) and photosystem II (PSII) stability/assembly 
factor HCF136 (spot L20), were identified in L. formosana. 
CAO, that converts chlorophyll a to chlorophyll b, regulates the 
stabilization of light-harvesting chlorophyll a/b proteins [26]. 
Villeth et al [27] reported that pathogen infection could decrease 
the accumulation of PSII stability/ assembly factor HCF136, 
which is important for the accurate assembly of PSII. In this 
study, AR remarkably decreased the abundance of CAO and PSII 
stability/ assembly factor HCF136 in L. formosana, but no light 
reaction-related proteins was depressed by AR in S. superba (Table 1 



and 2). These results suggest that photosynthesis apparatus of L. 
formosana is more sensitive to AR stress than S. superba. 

It has been reported that the expression of Calvin cycle enzymes 
were down-regulated in Arabidopsis under salinity stress [25] . Our 
previous study also reported that the reduction in photosynthesis 
was linked to Calvin cycle enzymes in AR- treated Arabidopsis [16]. 
Consistent with previous results, our proteomic data from this 
study confirmed that the abundances of Calvin cycle-related 
proteins including phosphoribulokinase (PPK, spot L28), Rubisco 
(spot L31, L34) and Rubisco activase (spot L27) were significandy 
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■= tj c decreased in L. fomwsana (Table 1). Rubisco, the CO2 fixing 

S o -S . . , .... 

E 5 % c enzyme in Calvin cycle, is the primary limiting factor of net 

S j£ xi m photosynthesis under environmental stress [16]. Rubisco activase 

~° To "5 I. promotes and maintains the catalytic activity of Rubisco [28]. The 

™ J S " decreased expression of Rubisco and Rubisco activase in AR- 

o h a p treated Z. fomwsana may disturb Calvin cycle activity, leading to 

the reduction in photosynthetic C0 2 assimilation and thus the 

_ inhibition in plant growth. Moreover, the results of Western blot 

§ -q g 3 analysis showed that the decrease in protein expression of Rubisco 

T3 large subunit was more obvious in L. formosana than in S. superba 

m c S. £ under AR stress (Fig. 6A and B). These results may explain why 
o -0 Q. S 

X 



X 



1 d 



s * g 

2 TO <- 
a- £ 



1 s 



ai £ 
2 x. 
> > 
S 53 



10 Q- 

i • 



c Xi aj 
- 5, ^ 



X 



AR-induced damage to photosynthesis was more serious in S. 
superba than in L. fmmosana (Fig. ID and E). Compared with L. 
formosana, less damage of AR to photosynthesis-related proteins 
« » J probably result from higher tolerance to AR in S. superba. 

£ S Besides photosynthesis-related proteins, AR also affected the 

I £ <2 abundances of energy production-related proteins in L. formosana 

% £ 8 and S. superba. It is well known that increased ATP production is 

a) z — required in response to abiotic stress in plants [29]. For example, 

u £ S the abundance of ATP synthase was considerably increased in 

J & ^ salt-stressed rice [30] and osmotic-stressed wheat [31]. Thus, the 

^ S. > increased abundances of ATP synthase subunits in L. formosana 

S| £ (spot L22, L23, L26) and S. superba (spot SI 1) in our study 

1 g "g demonstrate the prime role of ATP synthase in the adaptation of 

two tree species to AR stress. 



Material metabolism-related proteins 

In L. formosana, AR affected a series of protein abundances 
involved in several metabolism processes, including nitrogen 
metabolism (spot L2, L6, LI 1), starch and sugar metabolism 
u § (spot LI, L9), lipid metabolism (spot L7, L8), secondary 

Sj FJ ~ metabolism (spot L39, L46), vitamin metabolism (spot L3) and 

^ 2 ^ lignin biosynthesis (spot LI 2) (Table 1). However, only two 

o o j nitrogen metabolism-related proteins (spot SI, S3) and one starch 

!, |i biosynthesis-related protein (spot S2) were remarkably induced by 

a. AR in S. superba (Table 2). It is clear that AR disturbed more 

||| metabolism processes in L. formosana than in S. superba. 

~ S ^ AR stress can change free amino acid levels and disturbs N 

metabolism in plants [19]. Cysteine synthase (CS, spot LI 1), which 
is responsible for the terminal step of cysteine biosynthesis, is a 
» S g 5 critical enzyme involved in environmental stress response in plants 

jv o S [12]- The abundance of CS was decreased in L. formosana (Table 1), 

< i. o. o indicating that metabolic processes related to cysteine biosynthesis 

might be strongly depressed by AR. In addition, two glutamine 

^ » « S metabolism-related proteins, such as glutamate dehydrogenase 

h I — § I (GDH, spot L6) and glutamine synthetase (GS, spot S3), were 

identified in L. formosana and S. superba, respectively. Abiotic 

E "5 <l S stresses, such as salinity, drought and metal toxicity, can up- 

<s 2 £ c § regulate GDH gene expression and enhance GDH activity in 

2 aj 8~ plants [32]. GS catalyzes ATP-dependent incorporation of 

g H §■ jj a ammonium into glutamate and other reduced N compounds 

" H • 2 '5 [33]. It has been found that GS accumulation was also stimulated 

X °- i/i ™ o 

& 8 £ .g q. g by salt and drought, and this helped improve the tolerance of 

< f ^ plants to stresses [25]. In agreement with previous findings, our 
« S c " study found that the abundance of GDH in L. formosana, as well as 
-a sg it 'm "d g that of GS in 5. superba, were increased (Table 1), suggesting that 
15 « -S ™ ■§ c GDH and GS play critical roles in N metabolic acclimation of 
X> ~ ° § « — plants when exposed to AR. 

I ijj ~°- "S 3 £ Abiotic stress can also affect starch synthesis [34] . In this 
01 . « S " a 5, study, two starch biosynthesis-related proteins including glucose- 
5 ? d q " D ^ 1 -phosphate adenylyltransferase (GPAT, spot LI) in L. formosana 
" 2 5 U § o 1 5 and granule-bound starch synthase (GBSS, spot S2) in S. superba 
j= §" £■ I s 5 o were identified (Tables 1 and 2). Interestingly, the responses of 
GPAT and GBSS to AR were different in two broad-leaf species. 
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Signal transduction 
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Hormone response 
Unknown proteins 




I I L. formosana 
S. superba 



—i — 
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Number of protein spots 



Number of protein spots 



Figure 4. Number of protein spots significantly changed in AR-treated L. formosana and S. superba. (A) Protein spots increased in their 

abundances. (B) Protein spots decreased in their abundances. 

doi:10.1371/journal.pone.0102532.g004 



The abundance of GPAT was decreased in L. formosana. However, 
GBSS, the only enzyme implicated in amylose synthesis [35], 
was increased in AR-treated S. superba. Basically, this result was 
consistent with the observation where GBSS activity and starch 
content in rice was found to increase under cold stress [36]. It is 
well known that starch is required to synthesize sucrose which 
serves as a carbon and energy source for plant growth and stress 



response [35]. Thus we believe that the increased abundance of 
GBSS may contribute to higher AR-tolerance in S. superba through 
enhancing starch synthesis and energy production. It also should 
be noted that AR stress changed the abundances of lipid 
metabolism-related proteins including glycine-rich protein 17 
(GRP17, spot L7) and stearoyl-acyl-carrier protein desaturase-like 
protein (SACPDLP, spot L8), as well as secondary metabolism- 



Hormone response (3) 4.1% 

Others (12) 16.2% 

Post-translational 
modification (5) 6.8% 

Transcription (9) 12.2% 
Signal transduction (5) 6.8% 




Material metabolism (12) 
16.2% 



Photosynthesis and energy 
production (18) 24.3% 



Stress and defense (10) 
13.5% 



B 



Others (6) 17.6% 

Post-translational 
modification (2) 5.9% 

Transcription (4) 1 1 .8% 

Signal transduction (4) 11.8% 




Material metabolism (4) 1 1 .8% 



Photosynthesis and energy 
Production (10) 29.4% 



Stress and defense (4) 1 1 .8% 



Figure 5. Functional classification of AR-responsive proteins in L. formosana (A) and S. superba (B). The proportion of identities in each 

functional group was the sum of this identity accounting for all protein quantities. 

doi:10.1371/journal.pone.0102532.g005 
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L. formosana S. superba 
CK AR CK AR 



Rubsico LSU 




L. formosana S. superba 




L. formosana S. superba 



L. formosana S. superba 



Figure 6. Western blot analysis showing the expression of three protein spots. (A) Expression of rubulose-1,5-bisphoshate carboxylase 
large subunit (RuBisco LSU), ascorbate peroxidase (APX) and glutathione-S-transferase (GST) in L. formosana and S. superba seedlings after AR 
treatment. Relative expression level of RuBisco LSU (B), APX (C) and GST (D) were analyzed with the Quantity One software. /J-actin was used as the 
internal control. Means with different letters indicate significantly difference (p<0.05) with regard to AR treatments. 
doi:10.1371/journal.pone.0102532.g006 



related proteins including flavanone 3-hydroxylase (F3H, spot L39) 
and chalcone synthase (CHS, L46) in L. formosana. The abundance 
of SACPDLP was increased, which may contribute to lipid 
synthesis and membrane integrity in AR-treated L. formosana [25]. 
F3H and CHS are two key enzymes that catalyze the biosynthesis 
of flavonoids and chalcones, both of which play critical roles in 
enhancing secondary metabolism under environmental stress [37]. 
The increased abundances of these proteins imply that secondary 
metabolism pathway may be activated to cope with AR stress in 
AR-sensitive species, L. formosana.. 

Stress defense-related proteins 

The majority of environmental stresses generate a secondary 
oxidative stress in plants [13]. Oxidative stress occurs when there is 
a serious imbalance between ROS production and antioxidant 
defense [24]. Overproduction of ROS induced by a number of 
adverse environmental factors can attack proteins, lipids, and 
nucleic acids [38]. Under long-term heavy metal or AR stress, 
significant accumulation of ROS was observed in previous studies 
[7,19,39]. For instance, Kovacik et al [38] observed the increased 
ROS in four Tillandsia species under 2 U.M Cd + treatment over 30 
days. In this study, H 2 0 2 and 0 2 '" content, and thus oxidative 
stress, induced by AR was remarkably higher in L. formosana than 
in S. superba (Fig. 2C and D). 

To avoid oxidative damage, plants developed an antioxidant 
system consisting of antioxidative enzymes as well as non-enzymatic 



antioxidants [38] . In L. formosana, three antioxidant-related proteins, 
including APX (spot L40), GST (spot L41) and class III peroxidase 
ATP32 (spot L44), were identified by proteomic analysis. APX plays 
an important role in scavenging H 2 0 2 from cells [34]. GST is also 
an important enzyme that counteracts cellular damage induced by 
oxidative stress [12]. Enhanced activity of APX and GST has also 
regularly been detected in plants after exposure to salt, cold, heavy 
metal, and heat stresses [40]. Our recently published work found 
that both APX and GST were increased in their abundance in an 
AR-sensitive conifer tree species, Piuns massoniana [41]. Similarity, 
APX (spot L40) and GST (spot L41) with higher abundances and 
expression levels (Table 1 and Fig. 6) were also found in AR-treated 
L. formosana in this study, implying that antioxidant defense system 
was provoked by AR in L. formosana. However, the expression levels 
of APX and GST were not changed in S. superba seedlings under AR 
treatment (Fig. 6). A likely reason is that there may be other 
pathways that can remove excessive ROS in S. superba. Interestingly, 
thioredoxin peroxidase (TPx, spot SI 6), which appears to be a 
key enzyme in H 2 0 2 detoxification [42], was found to have an 
increased abundance in AR-treated S. superba. The increased 
TPx may function in resisting AR-induced oxidative damage by 
reducing ROS production in S. superba. In addition, endochitinase, 
a glycosyl hydrolase that catalyzes chitin degradation, plays an 
essential role in forming the fine cell-wall matrix that enhances the 
physical barrier against abiotic stresses [43]. Tapia et al [44] 
reported that endochitinase could be induced by heat, drought, and 
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salinity. In this study, the increased abundance of endochitinase 
(spot S 1 8) in S. superba (Table 2) may improve physical interactions 
at the plasma membrane-cell wall interface to cope with AR stress. 

Signal transduction-related proteins 

Signal transduction plays a crucial role in triggering a cascade of 
defense events [12]. Phytochrome is a chromoprotein that 
regulates the expressions of a large number of light-responsive 
genes and controls plant growth and development [45] . Recently, 
phytochrome has been found to modulate both biotic and abiotic 
stresses, such as salinity, drought, cold or herbivory [46]. Gross- 
talk between phytochrome-mediated light signals and some stress 
signaling pathways has been reported in diverse plants [47]. Thus 
it is possible that phytochrome is involved in the modulation of 
AR stress. A recent study found that increased abundance of 
phytochrome was needed to resist cold stress in cucumber [48]. 
Likewise, AR also increased the abundance of both phytochrome 
C (spot L48) and truncate phytochrome A2 protein (spot L49) in L. 
formosana but not in S. superba in this study (Table 1), which suggest 
a role of phytochrome signaling in response to AR in this sensitive 
species. In addition to the phytochrome signaling pathway, 
antioxidant enzymes have been found to be modulated by 
phytochromes under stress conditions [46]. In this study, the 
increased expression of antioxidant enzymes (APX and GST) was 
observed in L. formosana under AR treatment (Fig. 6). We suggest 
that the enhancement of phytochromes induced by AR modulates 
the antioxidant system in L. formosana seedlings. Further research 
need to widen our understanding of the role of phytochrome in 
AR-stressed plants. 

Calcium (Ca) also plays a crucial role in regulating plant 
defense responses to various environmental stimuli [49] . Recently, 
Kovacik et al [50] found that oxidative stress evoked by 
hexavalent chromium were evidently suppressed by Ca in 
Matricaria chamomilla using microscopic visualization method. Our 
previous study also reported that Ca addition dramatically 
alleviated the negative effects of AR on seed germination, seedling 
growth and photosynthesis [3]. Free Ca ion within cell is a second 
messenger for conveying internal and external signals by Ca 
sensors that subsequendy regulate diverse cellular processes in 
plants [17]. Calmodulin (CaM) and calcium-dependent protein 
kinase (CDPK), two major types of Ca sensor, play important roles 
in Ca signaling and further response to diverse stresses in plants 
[51]. Saijo et al [52] found that over-expression of OsCDPK gene 
enhanced tolerance to cold, salt and drought in transgenic rice. 
Moreover, 14-3-3 protein can also regulate proteins involved in 
stress response and activate CDPK signal transduction pathway in 
plants [13]. In S. superba, the abundances of CaM (spot S27), 
CDPK (spot S28) and 14-3-3 protein (spot SI 5) were increased 
after AR treatment, while these proteins were not identified in 
L. formosana. In agreement with the proteomic analysis results, 
the enhanced gene expression of CaM and CDPK by real-time 
quantitative PCR was more obvious in S. superba than in L. 
formosana under AR treatment (Fig. SI). These results indicate that 
AR activated different signaling transduction pathways in two tree 
species and that Ca sensor-dependent signaling pathway might 
play a critical role in enhancing AR tolerance in S. superba. 

AR affected transcription, protein synthesis and 
modification 

Transcription machinery plays an important role in abiotic 
stress adaptation [25]. Nine proteins (spots L52-L60) and four 
proteins (spots S23-S26) that involve in gene transcription were 
identified in response to AR stress in L. formosana and S. superba, 
respectively (Tables 1 and 2). Four of these proteins were maturase 



K (spots L57-59, S24). In plants, maturase K catalyzes intron 
RNA binding during reverse transcription and spicing and directiy 
affects gene expression at transcriptional level [53]. The changes in 
protein expression of maturase K are very complex in plants under 
abotic stress, depending on plant species and the type of stress 
[54]. Pandey et al [54] reported that maturase K was induced 
when plant suffered from high ROS pressure. On the contrary, 
maturase K protein was down-regulated in salt-treated maize [55]. 
In this study, the abundance of maturase K (spots L57-59) was 
increased in L. formosana, but was decreased in S. superba (Tables 1 
and 2). Based on our physiological data, AR induced higher ROS 
(H 2 0 2 and O2* ) production in L. formosana than that in S. superba 
(Fig. 2), which might be one reason to explain the increased 
abundance of maturase K in L. formosana. However, the decreased 
abundance of maturase K in S. superba indicates that its gene 
transcription had been affected by AR, thought no significant 
visible damage symptom emerged in S. superba leaves. 

Translational elongation factor Tu (EF-Tu) plays important 
roles in response to abiotic stresses including high and low 
temperatures, salinity, and water deficit [56]. Pandey et al [54] 
reported that the expression of EF-Tu was down-regulated in 
chickpea with the increased treatment period of water deficit, 
which was consistent with the decrease in EF-Tu abundance 
(spot LI 5) after AR stress in L. formosana observed in our study. 
This result indicates that AR has induced a lower protein synthesis 
in L. formosana. However, EF-Tu was not identified in S. superba, 
instead, proteasome subunit alpha type-5 isoform 1 (spot S4) and 
mitochondrial import receptor subunit TOM6 homolog isoform 1 
(spot S5) were found in AR-treated leaves (Table 2). Proteasome 
degrades the proteins with translational errors and proteins 
damaged by stress which can aggregate and become toxic to the 
cell [57]. Proteasome has been implicated in regulating numerous 
plant signaling and metabolic pathways under stress condition 
[57]. Our pervious study found that the expression of 20 S 
proteasome subunit was up-regulated in Arabidopsis after 32 h of 
AR treatment [16]. In this study, the abundance of proteasome 
subunit alpha type-5 isoform 1 was also increased in AR-treated S. 
superba, indicating that the control of protein degradation by the 
proteasome is likely to play an important role in enhancing AR- 
tolerance in S. superba. In addition, mitochondrial import receptor 
subunit TOM complex mediates the translocation and uptake of 
nuclear-encoded mitochondrial preproteins from the cytosol [58]. 
Increased abundance of mitochondrial import receptor subunit 
TOM6 by AR may accelerate the import of mitochondrial 
preproteins and enhance cellular metabolism and energy produc- 
tion in mitochondria, finally contributing to improve resistant to 
AR in S. superba. 

Hormone response-related proteins 

Plant hormones are not only important in plant growth and 
development, but also closely related to environmental stresses 
response [59]. It is known that increased ethylene biosynthesis is a 
general response of plants to stress conditions [60] . S-adenosyl- 
L-methionine (SAM), a major methyl donor in plants, is used as a 
substrate for ethylene biosynthesis [60]. SAM synthetase catalyzes 
SAM biosynthesis from L-methionine and ATP [61]. It has been 
found that both activity and expression of SAM synthetase were 
increased in tomato and rice under salt stress [62]. Accordingly, 
our present results indicate that AR increased the abundance of 
SAM synthetase (spot L5) in L. formosana (Table 1), which may 
further contribute to ethylene biosynthesis. Another ethylene- 
related protein, ethylene-responsive transcriptional coactivator 
(spot L42), which was also increased by AR stress in L. formosana, 
can positively control the expression of ethylene-responsive genes 
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in plants [63]. In addition, ABA inducible protein (spot L10), 
which is involved in ABA stimulus and abiotic stresses response, 
was increased in L. formosana when exposed to AR. It appeared 
that ethylene biosynthesis and ABA signaling pathways might be 
activated in AR-treated L. formosana. It is interesting to note that, in 
S. superba, no hormone response-related proteins was identified, 
partly due to the high resistance to AR stress in S. superba. 

Conclusion 

Using the approach of proteomic analysis, this study investigat- 
ed the differential responses to AR in two broad-leaf tree species, 
L. formosana and S. superba, an AR-sensitive species and an AR- 
tolerant species, respectively. After AR treatment, more proteins 
were significandy changed in their abundances in L. formosana than 
in S. superba. It should be noted that hormone response-related 
protein was only found in L. formosana. After AR treatment, 
signaling pathways, energy production and antioxidant system 
were activated in both L. formosana and S. superba. Due to higher 
AR-tolerance in S. superba, AR induced less damage to photosyn- 
thesis-related proteins in this species. Moreover, the proteins 
related to starch synthesis and translation were depressed in AR- 
treated L. formosana, but enhanced in S. superba, implying that these 
proteins may greatly contribute to enhance AR-tolerance in S. 
superba. The identification of novel AR-responsive proteins in this 
study provides not only new insights into AR stress responses, 
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